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Abstract 

We prove a double-exponential upper bound on the degree and on the 
complexity of constructing a Janet basis of a D-module. This generalizes 
a well known bound on the complexity of a Grobner basis of a module 
over the algebra of polynomials. We would like to emphasize that the 
obtained bound can not be immediately deduced from the commutative 
case. 



Introduction 

Let A be the Weyl algebra F[X\, . . . , X n , g^-, • ■ ■ , gjjH ( or the algebra of 
differential operators F{X\, . . . , X n )[-^j^, . . . , ^^H). Denote for brevity Di = 
oY - , 1 < i < n. Any A-module is called Z?-module. It is well known that 
an A-module which is a submodule of a free finitely generated ^4-module has a 
Janet basis. Historically, it was first introduced in [5]. In more recent times of 
developing computer algebra Janet bases were studied in [5], [13], [10]. Janet 
bases generalize Grobner bases which were widely elaborated in the algebra of 
polynomials (see e. g-[3]). For Grobner bases a double-exponential complexity 
bound was obtained in [15] , relying on pQ and which was made more precise 
(with a self-contained proof) in [3]. 

Surprisingly, no complexity bound on Janet bases was established so far; in 
the present paper we fill this gap and prove a double-exponential complexity 
bound. On the other hand, a double-exponential complexity lower bound on 
Grobner bases [T2], [H] provides by the same token a bound on Janet bases. 

There is a folklore opinion that the problem of constructing a Janet basis 
is easily reduced to the commutative case by considering the associated graded 
module, and, on the other hand, in the commutative case [6], [12], [4] the double- 
exponential upper bound is well known. But it turns out to be a fallacy! From 
a known system of generators of a D-module one can not obtain immediately 
any system of generators (even not necessarily a Grobner basis) of the associ- 
ated graded module. The main problem here is to construct such a system of 
generators of the graded module. It may have the elements of degrees {dl) 2 ° {n) , 
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see the notation below. Then, indeed, to the last system of generators of big 
degrees one can apply the result known in the commutative case and get the 
bound ((dlf W ) 2 uw = (dO 2 

. So new ideas specific to non-commutative 

case are needed. 

We are interested in the estimations for Janet bases of A-submodules of A 1 . 
The Janet basis depends on the choice of the linear order on the monomials (we 
define them also for I > I). In this paper we consider the most general linear 
orders on the monomials from A 1 . They satisfy conditions (a) and (b) from 
Section [1] and are called admissible. We prove the following result. 

THEOREM 1 For any admissible linear order on the monomials from A 1 any 
A-submodule I of A 1 generated by elements of degrees at most d ( with respect 
to the filtration in the corresponding algebra, see Section [7] and Section has 
a Janet basis with the degrees and the number of its elements less than 



We prove in detail this theorem for the case of the Weyl algebra A. The proof 
for the case of the algebra of differential operators is similar. It is sketched 
in Section [51 ^From Theorem [T] we get that the Hilbert function H(I,m), see 
Section Q] of the A-submodule from this theorem is stable for m > (dl) 2 ° 
and the absolute values of all coefficients of the Hilbert polynomial of I are 
bounded from above by (dl) 2 ^ ' , cf. e.g., [T3]. This fact follows directly from 
(fTT)|) . Lemma H2] from Appendix 1, Lemma [Hand Theorem [H We mention that 
in [7] the similar bound was shown on the leading coefficient of the Hilbert 
polynomial. 

Now we outline the plan for the proof of Theorem [T] The main tool in the 
proof is a homogenized Weyl algebra h A (or respectively, a homogenized alge- 
bra of differential operators h B) . It is introduced in Section [3J (respectively, 
Section [9|). The algebra h A (respectively h B) is generated over the ground 
field F by X , . . . , X n , Di, . . . , D n (respectively over the field F(X\, . . . , X n ) 
by Xq, Di, . . . , D n ). Here Xq is a new homogenizing variable. In the algebra A 
(respectively h B) relations (|T2|) Section [3] (respectively ([5fl|) Section [5]) hold for 
these generators in h A. 

We define the homogenization h I of the module I. It is a 'M-submodule of 
h A . The main problem is to estimate the degrees of a system of generators of I. 
These estimations are central in the paper. They are deduced from Theorem [2] 
Section [7] This theorem is devoted to the problem of solving systems of linear 
equations over the ring h A; we discuss it below in more detail. 

The system of generators of h I gives a system of generators of the graded 
gr(A)-module gr(I) corresponding to /. But gr(A) is a polynomial ring. Hence 
using Lemma 1121 Appendix 1 we get a double-exponential bound (dl) 20i ' on 
the stabilization of the Hilbert function of gr(7) and the absolute values of the 
coefficients of the Hilbert polynomial of gr(7). Therefore, the similar bound 
holds for the stabilization of the Hilbert functions of / and the coefficients of 
the Hilbert polynomial of I, see Section [2] 

But the Hilbert functions of the modules I and h I coincide, see Section [3] 
Hence the last bound holds also for the stabilization of the Hilbert functions of 
h I and the coefficients of the Hilbert polynomial of h I. In Section[5]we introduce 
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the linear order on the monomials from h A l induced by the initial linear order 
on the monomials from A 1 (the homogenizing variable Xo is the least possible 
in this ordering). Further, we define the Janet basis of I with respect to the 
induced linear order on the monomials. Such a basis can be obtained by the 
homogcnization of the elements of a Janet basis of / with respect to the initial 
linear order, see Lemma [3] 

Let Kdt( h I) be the monomial module (i.e., the module which has a system 
of generators consisting of monomials) generated by the greatest monomials of 
all the elements of the module h I, see Section |U Let C I, see Section HI be the 
module over the polynomial ring C A = F[X , . . . , X n , D\, . . . , D n ] generated by 
all the monomials from Hdt( h I) (they are considered now as elements of °A). 
Then the Hilbert functions of the modules h I and C I coincide. Thus, we have the 
same as above double-exponential estimation for the stabilization of the Hilbert 
functions of C I and the coefficients of the Hilbert polynomial of °I, Now using 
Lemma[T3]we get the estimation (d,l) 2 ° { ' on the monomial system of generators 
of C I, hence also of Kdt( h I). This gives the bound for the degrees of the elements 
of the Janet bases of I and hence also for the required Janet basis of /, and 
proves Theorem [T] 

The problem of solving systems of linear equations over the homogenized 
algebra is central in this paper, see Theorem O It is studied in Sections 0- 
[7J A similar problem over the Weyl algebra (without a homogenization) was 
considered in [7J . The principal idea is to try to extend the well known method 
due to G.Hermann [5] which was elaborated for the algebra of polynomials, to 
the homogenized Weyl algebra. There are two principal difficulties on this way. 
The first one is that in the method of G.Hermann the use of determinants is 
essential which one has to avoid dealing with non-commutative algebras. The 
second is that one needs a kind of the Noether normalization theorem in the 
situation under consideration. So it is necessary to choose the leading elements 
in the analog of the G.Hermann method with the least ordjsf , where Xo is a 
homogenizing variable, see Section [3J 

The obtained bound on the degree of a Janet basis implies a similar bound on 
the complexity of its constructing. Indeed, by Corollary [T] (it is formulated for 
the case of Weyl algebra but the analogous corollary holds for the case of algebra 
of differential operators) one can compute the linear space of all the elements 
z £ I of degrees bounded from above by (dl) 2 ° {n) . Further, by Theorem [1] the 
module Hdt(I), see Section[TJ is generated by all the elements Hdt(z) with z £ I 
of degrees bounded from above by (dl) 2 1 . Hence one can compute a system 
of generators of Hdt(I) and a Janet basis of I solving linear systems over F of 
size bounded from above (dl) 2 " 1 ' (just by the enumeration of all monomials 
of degrees at most (dl) 2 ° (n) which are possible generators of Hdt(/)). If one 
needs to construct the reduced Janet basis it is sufficient to apply additionally 
Remark Q] Section [4j 

For the sake of self-containedness in Appendix 1, see Lemma [121 we give a 
short proof of the double-exponential estimation for stabilization of the Hilbert 
function of a graded module over a homogeneous polynomial ring. A conversion 
of Lemma also holds, see Appendix 1 Lemma H31 It is essential for us. 
The proof of Lemma Q2] uses the classic description of the Hilbert function of 
a homogeneous ideal in F[Xq, . . . , X n ] via Macaulay constants b n +2, and 
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the constant bo introduced in [J. In Appendix 2 we give an independent and 
instructive proof of Proposition [T] which is similar to Lemma HUl In some sence 
Proposition [T] is even more strong than Lemma [13] since to apply it one does 
not need a bound for the stabilization of the Hilbert function. Of course, the 
reference to Proposition [T] can be used in place of Lemma fl3l in our paper. 

1 Definition of the Janet basis 

Let A — F[Xx, . . . , X n , D±, . . . , D n ], n > 1, be a Weyl algebra over a field F of 
zero-characteristic. So A is defined by the following relations 

X V X W = X W X V , D V D W = D W D V , D V X V —X V D V = 1, X V D W = D W X V , d/io. 

(1) 

By (|TJ) any element f £ A can be uniquely represented in the form 

/= E /< - , -v; ••••v; i >\ ■■■!>:;■ m 

*l.....«n,Jl,...Jn>0 

where all /ii,...,i„,ji,...,i n £ F and only a finite number of fi 1 ,...,i n ,j 1 ,...,j n are 
nonzero. Denote for brevity Z + ={z£/ : z > 0} to the set of all nonnegative 
integers and 

* V*l: ■ ■ • j*n): J = (jli ■ ■ • j jn)i fi.j — fi 1 ....,i n ,j 1 ....,j n 

X^X^...X^, I)- IH ...IP, . I V,,/, ; .V '/)■•. (3) 
\i\ =ii + ... + i n , i + j= (il +ji, ...,i„ +jn)- 

So i,j e Z" are multiindices. By definition the degree of / 

deg/ = deg Xl v .!>: n / = max{|i| + \ j\ : f lyJ ^ 0}. 

Let M be a left A- module given by its generators mi, . . . , m;, £ > 0, and relations 

E av, w m w , 1 <v < k. (4) 

l<u»<i 

where fc > and all a VjW £ A. We assume that dega VyW < d for all v, w. By ^ 
we have the exact sequence 

A k ^ A 1 ^ M -» (5) 

of left A-modules. Denote I = b(A k ) c A 1 . If I = 1 then J is a left ideal of A 
and M = A/I. In the general case I is generated by the elements 

(a B) i, . . . , a^) eA 1 , 1 < v < fc. 

For an integer m > put 

A m = {a : dega < m}, M m = ir(A l m ), I m = If] A l m . (6) 

So now A, M, I are filtered modules with nitrations A m , M m , I m , m > 0, 
respectively and the sequence of homomorphisms of vector spaces 

-> I m -> AL -» M m -> 
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induced by ([5]) is exact for every to > 0. The Hilbert function H(M,m) of the 
module M is defined by the equality 

H(M, to) = dim F M m , to > 0. 

Each element of A 1 can be uniquely represented as an F-linear combination of 
elements e v ^j = (0, . . . , 0, X l D 3 , 0, . . . , 0), herewith i,j G Z™ are multiindices, 
see (j3|), and the nonzero monomial X l D : > is at the position v, 1 < v < I. So 
every element / G A 1 can be represented in the form 

/ ^ ^ fv,i,j&v,i,j 5 fv,i,j G (7) 

The elements e^^j will be called monomials. 

Consider a linear order < on the set of all the monomials e Vl i : j or which is 
the same on the set of triples (v,i,j), 1 < v < I, i,j G Z™. If / ^ put 

o(/) = max{(v,i,j) : f v , itj ^ 0}, (8) 

see ©. Set 

o(0) = -oo < o(/) 

for every ^ / G A. Let us define the leading monomial of the element 
y£ f G A 1 by the formula 

Hdt^) — fv,i,j^v,i,ji 

where o(/) = j). PutHdt(O) = 0. Hence o(/-Hdt(/)) < o(/) if / ^0. For 
/ij/2 G A' if o(/i) < o(/ 2 ) we shall write /1 < / 2 . We shall require additionally 
that 

(a) for all multiindices for all 1 < v < I if i\ < i[, . . . ,i n < i' n and 
3i < j{, ■ ■ ■ , jn < in then (v,i, j) < (v,i',j r ). 

(b) for all multiindices i, j, for all 1 < v, v' < I if (v, i,j) < (v' , i', j') 
then (v,i + i",j + j")<(v',i' + i",j'+j"). 

Conditions (a) and (b) imply that for all /1 , / 2 G A 1 for every nonzero a G A 
if fi < Ji then a/i < a/ 2 , i.e., the considered linear order is compatible with 
the products. Any linear order on monomials e v jj satisfying (a) and (b) will 
be called admissible. 
Set 

Hdt(I) = J^HdtCf). 
fei 

So Hdt(J) is an ideal of A. By definition the family /1, . . . , f m of elements of / 
is a Janet basis of the module / if and only if 

1) Hdt(J) = AHdt(/i) + ...+AHdt(/ m ), i.e., the submodule of A 1 generated 
by Hdt(/i), . . . , Hdt(/ m ) coincides with Hdt(I). 

Further, the Janet basis /1, . . . , /,„ of / is reduced if and only if the following 
conditions hold. 

2) ft, . . . , f m does not contain a smaller Janet basis of J, 
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3) Hdt(A)>...>Hdt(/ m ). 

4) the coefficient from F of every monomial Hdt(/„), 1 < v < I, is 1. 

5) Let f a = J2 v .i,j fa,v,i,j&v,i,j be representation flS} for f a , 1 < a < to. 
Then for all 1 < a < (3 < m for all 1 < v < I and multiindices i, j the 
monomial f a , v ,i,j e v,i,j & Hdt(A/a \ {0}). 

Since the ring A is Noetherian for considered / there exists a Janet basis. Further 
the reduced Janet basis of I is uniquely defined. 

2 The graded module corresponding to a D— mo- 
dule 

Put A v = I v = M v = for v < and 

gr(A) = © m >oA m /A m _i, gr(7) = ® m > I m /I m -i, gr(M) = © m > M m /M m _i. 

The structure of the algebra on A induces the structure of a graded algeb- 
ra on gr(A). So we have gr(A) = F[Xi, . . . , X n , D\, . . . , D n ] is an algebra of 
polynomials with respect to the variables Xi, . . . , X n , D±, . . . , D n . Further, 
gr(/) and gr(M) are graded gr(A)-modules. From © we get the exact sequences 

-» Im/Im-l -» (4/4-i)' -» M m /M m _i -+ 0, TO > 0. (9) 

The Hilbert function of the module gr(M) is defined as follows 

ff(gr(M),m) = dim F M m /M m _ l7 m > 0. 

Obviously 

H(M,m)= #(g r ( M )> u ): i?(gr(M), m) = H(M, m) — H(M, m — 1). 

< 1! < m 

(10) 

for every to > 0. 

Denote for an arbitrary aGMby gr(a) £ gr(Af) the image of a in gr(M). 

LEMMA 1 Assume that b\, . . . , b s is a system of generators of I. Let Vi — 
deg hi, 1 < % < s. Suppose that for every to > 

I m = { X! Cvbv : Cv e ^' degc„ < to - 1 < i < s\. (11) 

1<D</J 

TTien gr(6i), . . . ,gr(6 s ) is a system of generators of the gr(A)-module gr(7). 
PROOF This is straightforward. 

So it is sufficient to construct a system of generators b\, . . . , b s oi I satisfying 

CD- 
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3 Homogenization of the Weyl algebra 

Let Xq be a new variable. Consider the algebra h A = F[Xo,Xi, . . . ,X n ,Di, 
. . . , D n ] given by the relations 

X V X W = X W X V , D V D W = D W D V , for all v,w, ,..„, 
D V X V - X V D V = Xg, 1 < v < n, X V D W = D W X V for all v w. 1 ' 

The algebra h A is Noetherian similarly to the Weyl algebra A. By (jTSJ) an 
element / £ h A can be uniquely represented in the form 

/= E / - , K ••••v; n\ ---IK ■ (13) 

io,ii,—,in,ji,—Jn>0 

where all fi ,...,i„,j-,_,...,j n £ F and only a finite number of fi ,...,i n ,j 1 ,...,j„ are 
nonzero. Let i,j be multiindices, see ([3]). Denote for brevity 

i — (il, . . . ,i n ), j = (Jx, . . . ,j n ), fi ,i,j = fi ,...,i n ,j 1 ,...,j n (1A.) 

/ v ,.v: .v^/>'. 1 j 

By definition the degrees of / 

deg/ = deg Xo . „ tXn ,D x ,.. .,D n f = max{i + \i\ + \ j\ : /i ,i,j 7^ 0}, 
deg Dl) ... jlJn / = max{U| : /i ,i,j + 0}, 
degr, Q / = max{j Q : / io;i j ^ 0}, 1 < a < n 
degx / = max{i Q : f io>i j ^ 0}, 1 < a < n 

Set ordO = ord Xfl = +00. If ^ f e h A then put 

ord / = ord Xo / = fi if and only if / £ Xg( h A) \ Xg +X ( h A), fi>0. (15) 

For every z — (zi, . . . ,z{) £ h A l put 

ordz = minjordzi}, degz = maxjdegZi}. 

l<i<l l<i<l 

Similarly one defines ord b and deg b for an arbitrary (fc x Z)-matrix b with 
coefficients from h A. More precisely, one consider here & as a vector with kl 
entries. 

The element / £ h A is homogeneous if and only if fi 0l ij ^ implies «o + \i\ + 
\j\ = deg/, i.e., if and only if / is a sum of monomials of the same degree deg /. 
The homogeneous degree of a nonzero homogeneous element / is its degree. 
The homogeneous degree of is not defined (0 belongs to all the homogeneous 
components of h A, see below). 

The m-th homogeneous component of h A is the i^-linear space 

( A) m = { z € A : z is homogeneous & deg z = m or z = } 

for every integer m. Now h A is a graded ring with respect to the homogeneous 
degree. By definition the ring h A is a homogenization of the Weyl algebra A. 

We shall consider the category of finitely generated graded modules G over 
the ring A. Such a module G = ® m >m G m is a direct sum of its homogeneous 
components G m , where m, tuq. are integers. Every G rn is a finite dimensional 
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.F-linear space and ( h A) p G m C G p + m for all integers p, m. If G and G' are 
two finitely generated graded M-modules then ip : G — > G' is a morphism (of 
degree 0) of the graded modules if and only if ip is a morphism of 'Vl-modulcs 
and (p(G m ) C G' m for every integer m. 

The element z G A (respectively z G A) is called to be the term if and 
only if z — Xz± ■ . . . ■ z v for some ^ A G -F, integer v > and z m G 
{X , . . . ,X n ,D 1 , . . . , D n ] (respectively z w G {X 1} . . . , X n , D x , . . . , D n }), 1 < 
w < v. 

Let z = J2j z j G ^4 be an arbitrary element of the Weyl algebra A represented 
as a sum of terms Zj and deg z = maxj deg Zj. One can take here, for example, 
representation (|3|) for z. Then we define the homogenization h z G h A by the 
formula 

h _ ■? vdegz-degZj 

z — ^ Zj^o 

3 

By (fTJ) , (fT^I the right part of the last equality does not depend on the chosen 
representation of z as a sum of terms. Hence h z is defined correctly. If z G h A 
then "z £ A is obtained by substituting Xq = 1 in z. Hence for every z G A we 
have < " l z = z, and for every z G ''A the element z — ,m zX^, where /i = ordz. 
For an element z = (zi, . . . , z\) G A 1 put degz = maxi<i<;{deg z{\ and 

,l z = ('^X^ 2 ^ 21 , . . . ^X^*-^* 1 ) G h A l . 

Similarly one defines deg a and the homogenization h a = (a v . w )i<v<k. i<w<i for 
an arbitrary k x /-matrix a with coefficients from A. More precisely, one consider 
here a as a vector with kl entries. Hence if b = (b VyW )i< v <k,i<tu<i — a then 
b v ,w = h o,v,wX^ csa ■^s^-"' f or a ii v,w. 

The m-th homogeneous component of h A l is 

( h A l ) m = { h z : z G A 1 & degz = m or z = } 

For an i^-linear subspace X C A 1 put ,l X to be the least linear subspace of h A l 
containing the set { h z : z G X}. If X is a (finitely generated) A-submodulo of 
A 1 then h X is a (finitely generated) graded submodule of h A l . The graduation 
on h X is induced by the one of h A . 

For an element z = (z%, . . . , zi) G M.' put a z = ( a zi, . . . , a z;) G A 1 . For a 
subset X C put Q X = { a z : z G X} C A'. If X is a F- linear space then a X 
is also a F-linear space. If X is a finitely generated graded submodule of h A l 
then °X is finitely generated submodule of A 1 . 

Now h I is a graded submodule of h A\ Further, "*/ = /. Let ( h I) m be the 
m-th homogeneous component of /. Then 

Vm) = ffio< J <™(^ 1 rn>0, (16) 
\{ h I) m ) = Im, m>0. (17) 

and (jTTJ) induces the isomorphism l : ( h I) m — > I m . Set "M = h A l / h I. Hence 
h M is a graded ^-module and we have the exact sequence 

-»■ h I -> 'M' -> ' l Af 0. (18) 

The m-th homogeneous component ( h M) m of 

( h M) m = ( h A l ) m /( h I) m - A l JI m . (19) 
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by the isomorphism l. Wc have the exact sequences 



-» ( h I) 



m 




( h M) 



rn 



o, m > o. 



(20) 



By definition the Hilbert function of the module M is 

H( h M, m) = dim F ( h M) m , m > 0. 

By (|19|) we have H(M,m) = H( h M,m) for every m > 0, i.e., the Hilbert 
functions of M and M coincide. 

LEMMA 2 Let b%, . . . ,b s be a system of homogeneous generators of the h A- 
module I. Then 



is a system of generators of gr(A) -module gr(/). 

PROOF By a {{ h I) m ) = I m . Now the required assertion follows from 
Lemma [TJ The lemma is proved. 

4 The Janet bases of a module and of its ho- 
mogenization 

Each element of A can be uniquely represented as an F-linear combination of 
elements e v<io>ilj = (0, . . . , 0, X l "X l D^ 0, . . . , 0), herewith < i G Z, i, j G 
are multiindices, see ((3]), and the nonzero monomial X 1 qX 1 D : > is at the position 
v, 1 < v < I- So every element / G can be represented in the form 



and only a finite number of f Vl i 0l ij are nonzero. The elements e v ,i ,i,j will be 
called monomials. 

Let us replace everywhere in Section [T] after the definition of the Hilbert 
function the ring A, the monomials e v ,i,ji the multiindices i, i' , i" , triples 
(v,i,j), (v,i',j'), the module / and so on by the ring "A, monomials e Vt i 0t % t ji 
the pairs (i' ,i'), (ig,i") (they are used without parentheses), quadruples 

(v,io,i,j), (v, i' 0) i\ j'), the homogenization h I and so on respectively. Thus, we 
get the definitions of o(f), Hdt(/) for / G A , new conditions (a) and (b) which 
define admissible linear order on the monomials of A , new conditions l)-5), 
the definitions of the Janet basis and reduced Janet basis of h I. For example, 
the new conditions (a) and (b) are 

(a) for all indices io>*o> au multiindices for all 1 < v < I if «o < *oj 
k <i[,---,in< i' n an d ji < j{, ■ ■ ■ ,j n < j'„ then (v,i Q ,i,j) < (v,i' , i', j'). 

(b) for all indices io,i' ,iQ, all multiindices for all 1 < v,v' < I 
if {v,i ,i,j) < (v',i' ,i',j') then (v,i + i'o,i + i",j+j") < (v',i' + i'o,i' + 



g r( < % 1 ),... s gr('% a )6gr(A) 




G F. 



(21) 



i",f + j") 
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The Janet basis of I is homogeneous if and only if it consists of homogeneous 
elements from h A . 

Let < be an admissible linear order on the monomials from A 1 , or which is the 
same, on the triples (v, i,j), see Section[T] So < satisfies conditions (a) and (b). 
Let us define the linear order on the monomials e v ,i ,i,j or, which is the same, 
on the quadruples (v,io,i,j). This linear order is induced by < on the triples 
(v,i,j) and will be denoted again by <. Namely, for two quadruples (v,io,i,j) 
and (v',i' ,i',f) put (v,i ,i,j) < (V, i' , i',f) if and only if (v,i,j) < (v',i',j'), 
or (t>, i, j) — (v' , z', j') but io < i' . Notice that this induced linear order satisfies 
conditions (a) and (b) (in the new sense). 

REMARK 1 If /i , . . . , f m is a Janet basis of I ( respectively homogeneous 
Janet basis of I ) satisfying 1)~4) then there are the unique c a ^ 6 A (respec- 
tively c at p £ A), 1 < a < (3 < ra, such that 

fa + c a,/3f/3, l<a<m, 

a < /3 < m 

is a reduced Janet basis of I ( respectively reduced homogeneous Janet basis of 
h I), cf 0/. 

LEMMA 3 Let fx, ... , f m be a (reduced) Janet basis of I with respect to the 
linear order <. Then h f 1 , . . . , h f m is a (reduced) homogeneous Janet basis of the 
module h I with respect to the induced linear order < . Conversely, let g\ , . . . , g m 
be a (reduced) homogeneous Janet basis of the module I with respect to the 
induced linear order <. Then a g 1 ,... , a g m is a (reduced) Janet basis of I with 
respect to the linear order <. 

PROOF This follows immediately from the definitions. 

Let / £ A and the module h I be as above. Then there is the unique element 
g e h A l such that 

.9 = ^ ] 9v,i ,i,j e v ,ia,i,j> 9v,ia,i,j 6 ^1 

f — g E h I and if g v ,i ,i,j ^ then e Vi i 0: ij ^ Hdt( h I). The element g is called the 
normal form of / with respect to the module I. We shall denote g — nf( h I, /). 
Obviously nf( /l J, {^A l ) m ) C { h A l ) m is a linear subspace. 

Let C A = F[Xq, . . . , X n , Di, . . . , D n ) be the polynomial ring in the variables 
Xq, . . . , X n , Di,..., D n . Each monomial e Vi i .ij can be considered also as an 
element of c A l . Denote by C I C c A l the graded submodule of c A l generated by 
all the monomials e V: i a ^j such that there is 7^ / 6 I with o(f) = (v,io,i, j). 
The Hilbert functions 

H( C I, m) = dimi?{(zi, . . . , zi) S C I : Vz ( degz^ = m or z,-=0)}, 

H(°A l /%m) = ( m ^ 2n ) - H(%m). 

Let us replace in the definition of the normal form above A, h I by C A, C I respec- 
tively. Thus, for / G c A l we get the definition of the normal form nf( c J, /) G c A l , 
cf. 0]. Obviously, nf( c 7, ( c A l ) m ) C { c A l ) m is a linear subspace. Since the ideals C I 
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and Kdt( h L) are generated by the same monomials we have dimnf ( C I, ( c A l ) m ) = 
dimnf(' 1 /, ( h A) m ). Hence the Hilbert functions 

H( h A l / h I, m) = H( c A l / c I, m), H( h I, m) = H{ C I, m), m > 0, 

coincide. Therefore, see Section [31 

H(I,m) = H( c I,m), m > (22) 

5 Bound on the kernel of a matrix over the ho- 
mogenized Weyl algebra 

LEMMA 4 Let k = I — 1 and I > 1 be integers. Let b — (£'i,j)i<i<fe, i<j<i be a 
matrix where bij G A are homogeneous elements for all Let degbij < d, 
d > 1, for all Assume that there are integers dj > 0, 1 < i < k, and d[ > 0, 
1 < j < I, such that 

degbij =di- d'j (23) 

for all nonzero bij, and additionally mini<j<i{dj} — (hence di < d, d'^ < d 
for all i,j), d > 1. Then there are homogeneous elements Zi,...,zi € A such 
that (zi, . . . ,zi) ^ (0, ... ,0), 

y^b i>j Zj = Q, 1 < * < Z — 1, (24) 
i<j<i 

all nonzero bijZj have the same degree depending only on i and 

deg zj < (2n + 3)ld, l<j<l. (25) 

Besides that, if all bij do not depend on X n (i.e., they can be represented as 
sums of monomials which do not contain X n ) then one can choose also z\,...,zi 
satisfying additionally the same property. Finally, dividing by an appropriate 
power of Xq one can assume without loss of generality that minjordzi : 1 < 
i < 1} = 0. 

PROOF We shall assume without loss of generality that I > 2. At first 
suppose that that degbij = deg b for all nonzero bij. Consider the linear 
mapping 

(hA\l , (hA\l-l 

(zi,...,z t ) h- (T,i<j<iKjZj) . 

\ — J — / Ktsl-1 



(26) 



l[ m - d l gb + 2n )>(l-l)( m + 2n ) (27) 



If 

jj'm — c 

2n J " v " ~'\ 'li> 

then the kernel of (|26p is nonzero. But l[27p holds if 

1 + ^ W 1+ deg5 V,,fi+ de l b W-L. 

77i + 2n — deg o / \ m + 2n — 1 — deg 6 / V m — deg b J I — 1 

(28) 

Further, (gSJ) is true if (1 + deg6/(m - deg6)) 2n < 1/(1 - 1). The last inequality 
follows from m > (2n + 1) degb/ log(Z/(Z — 1)). Hence also from m > (2n + 
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l)Zdeg&. Notice that (2n + 2)ld > 1 + (2n + l)ldegb. Thus, the existence of 
Zi, . . . , zi is proved, and even more all nonzero bijZj have the same degree which 
does not depend on i,j. Notice that in the considered case we prove a more 
strong inequality deg Zj < (2n + 2)ld for all 1 < j < I, 

Suppose that ax, . . . , a; do not depend on X n . We represent Zi = J2j z i,j^-rn 
1 < i < I, where all Zij do not on X n . Let a = max.;{dcg Xii z{\. Obviously in 
this case one can replace (z\, . . . , z{) by {z\ >a , . . . , z^ a ). 

Let us return to general case of arbitrary deg bij. We shall reduce it to 
the considered one. Namely, multiplying the i-th equation of system (|24p to 
^max,{d,} we sna u suppose without loss of generality that all dk are equal. 

Let us substitute ZjX 3 for zj in ([24| . Now the degrees of all the nonzero 
coefficients of the obtained system coincide. Thus, we get the required reduction 
and estimation (|25[) . The lemma is proved. 

REMARK 2 Lemma^ remains true if one replaces in its statement condition 
^by 

^2 z jh,j =0, 1 < i < I - 1, (29) 

1<7"<Z 

The proof is similar. 

REMARK 3 Let the elements bij be from Lemma [7J Notice that there are 
integers 5[ > 0, 1 < i < k, and Sj > 0, 1 < j < I, such that 

deg =Sj -S'i 

for all nonzero bij, and mini<i<fc{(5,'} = 0. Namely, 5[ = —di + maxi<i<k{di}, 
Sj = -d'j + maxi<i<fc{<ii}. 

6 Transforming a matrix with coefficients from 
h A to the trapezoidal form 

Let b be the matrix from Lemma U] but now k, I are arbitrary. Hence (|23p 
holds. Let b = (bi, ...,&;) where b\, . . . , h £ h A k be the columns of the matrix 
b (notice that in Lemma [1] and Lemma [2] bi are rows of size I; so now we change 
the notation). By definition bi, . . . , bi are linearly independent over h A from the 
right (or just linearly independent if it will not lead to an ambiguity) if and only 
if for all zx, ■ ■ ■ , Zi S A the equality bxzx+- ■ -+hzi = implies z\ = . . . = Zi = 0. 
By (|23|) in this definition one can consider only homogeneous zx , ■ • ■ , zi . For an 
arbitrary family &i, . . . , bi from Lemma |4] (with arbitrary fc, I) one can choose a 
maximal linearly independent from the right subfamily bi 1 , . . . , bi r of b\ , . . . , bi . 
It turns out that r does not depend on the choice of a subfamily. More precisely, 
we have the following lemma. 

LEMMA 5 Let Cj = X)i<i<i^ Zi J' ^ — 3 — r ^> w ^ ere z i,j S A are homo- 
geneous elements. Suppose that there are integers d'L 1 < j < rx, such that 
for all i,j the degree degZjj = d[ — d'j. Assume that Cj, 1 < j < rx, are lin- 
early independent over A from the right. Then r\ < r, and if rx < r there are 
cv-l+i, . . . , c r 6 {bi ± , . . . ,bi r } such that Cj, 1 < j < r, are linearly independent 
over h A from the right. 
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PROOF The proof is similar to the case of vector spaces over a field and we 
leave it to the reader. 

We denote r = rankr{&i, . . . ,bi} and call it the rank from the right of b-y, ■ ■ . , 6;. 
In the similar way one can define rank from the left of b\ , . . . , bi . Denote it by 
rankl{&i, . . . , b{\. It is not difficult to construct examples when rankr{6i, . . . , bi] 
^ rankl{&i, . . . , b{\. The aim of this section is to prove the following result. 

LEMMA 6 Let b be the matrix with homogeneous coefficient from h A satisfying 
i23\) . see above. Suppose thatdegbij < d for alii, j. Assume thatk > I > 1. Let 
li = rankr{&i, . . . , b{\ and b\, . . . , bi t be linearly independent. Hence < l± < I. 
Then there is a matrix {zj, r )i<j,r<h with homogeneous entries Zj, r € h A and a 
square permutation matrix a of size k satisfying the following properties. 

(i) All the nonzero elements bijZj_ r for 1 < j < I have the same degree 
depending only oni,r and 

degz j>r <{2n + 3)ld. (30) 



(ii) Set the matrix e = (eij)i<,-<fe, i<j<h = &bz. Then the matrix 

e' 
e" 

where e' = diag(e' 1 1 , . . . ,e[ . ) is a diagonal matrix with l\ columns and 
each e'j j, 1 <j < l±, is nonzero. 

(Hi) orde^j > orde^ for all 1 <i < k, I < j <l\. 

Besides that, if all (and hence all bi,j) do not depend on X n (i.e., they 
can be represented as sums of monomials which do not contain X n ) then one 
can choose also Zj, r satisfying additionally the same property. Finally, dividing 
by an appropriate power of Xq one can assume without loss of generality that 
minlordzj,,. : 1 < J < ii} = for every 1 < r < l\. 

PROOF At first we shall show how to construct z and e such that (ii) and 
(iii) hold. We shall use a kind of Gauss elimination and Lemma [U Namely, we 
transform the matrix e. At the beginning we put 

e = (ei, . . .,e h ) = . . .,b h ). 

We shall perform some M-linear transformations of columns and permutations 
of rows of the matrix e and replace each time e by the obtained matrix. These 
transformation do not change the rank from the right of the family of columns 
of e. At the end we get a matrix e satisfying the required properties (ii), (iii). 

We have rankr(e) = l\. If l\ = 0, i.e, e is an empty matrix, then this is the 
end of the construction: z' is an empty matrix. Suppose that l\ > 0. Let us 
choose indices 1 < io < k, 1 < jo < l\ such that ordei ojo = mini^xjjlordej}. 
Permuting rows and columns of e we shall assume without loss of generality 
that {i ,j ) = (1, 1). 

By Lemma H] we get elements w^\,w^i E h A of degrees at most (2n + 3)2d 
such that ei t iWi,i — ei^w^i, 1 < i < h, and ordu^i = for every 1 < i < l\. 
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Set w' = (—Wi t 2, ■ ■ ■ , —W\,i x ), and w" = diag(w2,2, ■ ■ ■ , Wh,h) t° ^ e the diagonal 
matrix. Put 

w= {o, Z") 

to be the square matrix with l\ rows. We replace e by ew. Now 

( E2,i, E 2 .2 ) ' 
where £^2,2 has l\ — 1 columns and 

min {ov&bj} = ordei.i = min {orde.,} (31) 
i<j<h " ' i<j<li 

(for the new matrix e). 

Let us apply recursively the described construction to the matrix £2,2 in 
place of e. So using only linear transformations of columns with indices 2, . . . , l\ 
and permutation of rows with indices 2, . . . , k we transform e to the form 

/ e M; 
erer = S 2 x , 

V eL 

where a is a permutation matrix and r' is a square matrix with l\ — 1 rows (it 
transforms £2,2)) the matrix E' 22 = diag(e2.2, ■ ■ • , is a diagonal matrix 

with l\ — 1 > columns, and all the elements e2,2, ■ ■ ■ >&h,h 6 ft A are nonzero. 
We shall assume without loss of generality that a = 1 is the identity matrix. We 
replace e by er. Conditions (ii) and (hi) hold for the obtained e and, more than 
that, by (iii) applied recursively for (E 2: 2, E' 22l E 22 ) (in place of (e, e', e")), and 
pip the same equalities are satisfied for the new obtained matrix e. 

Let E'21 = (e 2 ,i, . . . , e/j^.i)* where t denotes transposition. By Lemma H 
there are nonzero elements 1*1,1, ■ ■ • , v h,i € A of degrees at most 

(2n + 3)(max{dege M : l<i<h} + l)h (32) 

such that e%jy\ i = e^jt^i and min{ord vi,!, ordwi^} = for all 1 < i < l\ — 1. 
Let v' — (—1)2,1, ■ ■ ■ j ~ v h,iY and v" be the identity matrix of size l\ — 1. Put 

f «i,i> ^ 

Let us replace e by ew. Put z — wtv, where the matrix z has l\ columns. Recall 
that without loss of generality a — 1 is the identity permutation. We have 
e = (61, ... , bi^z. These Gauss elimination transformations of e do not change 
the rank from the right of the family of columns of e. It can be easily proved 
using the recursion on I, cf. Lemma[5]below. Now the matrix e satisfies required 
conditions (ii), (iii) and a = 1. 

Let us change the notation. Denote the obtained matrix z by z' . Let z' = 
(z^, . . . , z 1 ^) where z'j is the j'-th column of z' . Our aim now is to prove the 
existence of the matrix z satisfying (i)-(iii). By Lemma |4] for every 1 < r < l\ 
there are homogeneous elements Zj jr € A, 1 < j ' < I, such that (zi t r, . . . , z^ r ) ^ 
(0,...,0), 

/J bi,jZj tr = for every I <i <h, i ^= r, (33) 

i<j<h 
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and estimations for degrees (J3UJ) hold. Put the matrix z — (zj, r )i<j,r<li- 

Let 

z = (zi, . . . , zi t ) where Zj is the j-th column of z. Hence Zj — (zi jr , . . . , z/. r )'. 
LEMMA 7 For every 1 < r < l\ we have 

Kj z i,r ± 0. (34) 

1<7<*1 

Further, for every 1 < r < l\ there are nonzero homogeneous elements g' r ,g r £ 
h A such that z' r g' r = z r g r . 

PROOF Consider the matrix (z', z r ) with l\ rows and l\ + 1 columns. By 
Lemma|3]there are homogeneous elements hi, ... , hi 1+ i 6 M (they depend on r) 
such that (hi,..., hi 1+ i) ^ (0, . . . , 0) and the following property holds. Denote 
h=(hi,...,h h+ if, h'= (hi,...,h h y. Then 

z'h! + z r h h+ i = (35) 

(we don't need at present any estimation on degrees from Lemma 21 only the 
existence of h). Denote by b" the submatrix consisting of the first rows of the 
matrix (bi, . . . , bi t ). Multiplying (|3"5"]) to b" from the left we get 

b"z'ti + b"z r h h+1 = 0. (36) 

But b" z' is a diagonal matrix with nonzero elements on the diagonal, see (ii) 
(for z' in place of z). Hence by ([3"3"|l and (|36p hj — for every j ^ r. Now 
h ^ (0,...,0)* implies h r ^ and h h+1 ^ 0. Therefore, (011) holds. Put 
g' r = h r and y r = hi 1+ \. We have z^.g^. = z r g r by ([55]). The lemma is proved. 

Let us return to the proof of Lemma [HI Now (i)-(iii) are satisfied by Lemma [7] 
The last assertions of Lemma E] are proved similarly to the ones of Lemma HI 
Lemma |6] is proved. 



7 An algorithm for solving linear systems with 
coefficients from h A. 

Let u = (tii, ■ ■ ■ , u iY S A . Let all nonzero be homogeneous elements of the 
degree — d'j + p for an integer p. Suppose that — d'j + p < a" for an integer d! > 1. 
Let & = (bij)i<i<k. i<j<i be the matrix with fc rows and I columns from the 
statement of Lemma [S] (but now k and I are arbitrary). So degfe^j = di — dj < d 
for all i, j. Let Z = (Zi, . . . , Zk) be unknowns. Consider the linear system 



y t Zikj =uj, l<j<l, (37) 

Ki<k 



or, which is the same, 
Denote 



Zb- 



ordu= min {orditi}. (38) 

l<i<k 

The similar notations will be used for other vectors and matrices. In this section 
we shall describe an algorithm for solving linear systems over h A and prove the 
following theorem. 
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THEOREM 2 Suppose that system |ff7p has a solution over A. One can 
represent the set of all solutions of Jff7] j over h A in the form 

J + z*, 

where J C A is a h A-submodule of all the solutions of the homogeneous system 
corresponding to {3ty (i.e., system f39j ) with all Uj = 0) and z* is a particular 
solution of Ij37\ ). Moreover, the following assertions hold. 

(A) One can choose z* such that ordz* > ordu — v , where v > is an integer 
bounded from above by (dl) 2 ° ' (and depends only on d and I ). The degree 
deg z* is bounded from above by d'(dl) 2 ( . 

(B) There exists a system of generators of J of degrees bounded from above by 
(dl) 2 ( . The number of elements of this system of generators is bounded 
from above by k(dl) 2 ° ( ' . 

Besides that, if all bij and Uj do not depend on X n (i.e., they can be represented 
as sums of monomials which do not contain X n ) then z* and all the generators 
of the module J also satisfy this property. 



PROOF Let l\ = rankr(&i, . . . , b{). Permuting equations of (j3"T|) we shall 
assume without loss of generality that (b%, . . . , b^) are linearly independent from 
the right over "A. Let a, z, e, e', e" be the matrices from Lemma[U Similarly to 
the proof of Lemma [6] we shall assume without loss of generality that a = 1 . 
Denote by b' the submatrix of b consisting of the first l\ columns of b, i.e., 
b' = (&i, ... ,bi x ). By Lemma H] there are nonzero elements qi t x, . . . , qi 1 j 1 of 
degrees at most such that e\^q\ t i = e.i t iq\i t i and minjordqi i, ordq^} = 
for all 2 < i < l\. Set q = diag(qi.i, . . . , q^.^) to be the diagonal matrix. Let 
vq = orde^igi,!. Then by Lemma [6] (iii) ord(6'zq) > v$. Let Xq°5 = b' zq. 
Then 8 is a matrix with coefficients from A and 



5 = 



5' 
5" 



where S' = diag(5i l i, . . . , Si lt i ± ) is a diagonal matrix with homogeneous coef- 
ficients from h A and all the elements on the diagonal are nonzero and equal, 
i.e., 5jj = <5i.i for every 1 < j < 1%. Besides that, ord^i = 0. Fur- 
ther, 8" — {5i.j)i 1 +i<i<k.i<j<i 1 - We have ord(uzq) > vq, since, otherwise, 
system ([57]) does not have a solution. Obviously ordu < ord(uzq). Denote 
u 1 = (uq, . . . ,u[Y — X^ v °uzq E A . Hence ordit' > ord(w) — Vq. Consider the 
linear system 

ZS = u'. (39) 

LEMMA 8 Suppose that system has a solution over A. Then linear 

system i39\) is equivalent to j37\ ), i.e., the sets of solutions of systems &39\) and 
over h A coincide. 

PROOF The system Zb 1 z — uz is equivalent to |37|) by Lemma El System 
(|3U) is equivalent to Zb' z = uz since the ring h A does not have zero-divisors. 
The lemma is proved. 
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REMARK 4 Since rankr(6i, . . . , b{) = 1% and by Lemma\fi\for every l\ + 1 < 
j < I there are homogeneous Zjj, zij, . . . , z\ x j e A such that Zjj ^ and 
bj z j,j + Si<r<2i b r z r .j = and all degZjj, degz rj - are bounded from above by 
(2n + 3)(Zi + ljd. Put u'j — UjZjj +Xa<r<fi u r z r,j, h + 1 < j < I- Then system 
([371 ) has a solution if and only if system \39\) has a solution and u'j — for all 
l\ + 1 < j < I. This follows from Lemma\$\ and Lemma\^[ But in what follows 
for our aims it is sufficient to use only Lemma\8[ 

REMARK 5 Assume that deg^ bij < for all i.e., the elements of the 
matrix b do not depend on X n . Then by Lemma^ and the described construction 
all the elements of the matrices b, z, q, S, 5', 6" also do not depend on X n . 

By Lemma and Remark [5] for every l± + 1 < j ' < k there are homogeneous 
elements 9j,j,gj,i S A, 1 < i < h, such that 

9jj$j,i = 9j,ih,u 1 < i < h, 

all the degrees deg gjj, deggr^i, 1 < i < l\, are bounded from above by 

(2n + 3)(Zi + l)(max{deg^, i : 1 < i < k} + 1) 

and mhiKi^ {ordtfyj, ord gj,i} = 0. Hence ord gj j — for every li + 1 < j < k 
since ord^i^ = 0. 

Denote h = 6i,ifl'j 1 +i.,j 1 +i<7ii+2,ii+2 • ■ ■ fffc.fc. So h e h A is a nonzero homoge- 
neous element and ord/i = 0. Set e = degh. We need an analog of the Noether 
normalization theorem from commutative algebra, cf. also Lemma 3.1 [7]. 

LEMMA 9 There is a linear automorphism of the algebra A 
a: h A^ h A, a(Xi) = •:'>, ,, A, • a 2 ,i, j D j ), 

l<j<n 

Ci(Di) = ^ ( a 3,i,j X j + a i,i,3 D i)i a ( X o) = X , 1 < i < U, 
l<J<n 

such that all ot s ,i,j S F, deg^ a(h) — e. If deg Xn h = then one can choose 
additionally a(X n ) = X n , all cti in< j — for 1 < j < n — 1 and ct3,n.j = for 
1 < j < n. 

PROOF Recall that ord h — 0. Hence at first it is not difficult to construct 
a linear automorphism (3 such that (3{Xq) = Xq, 

(3{X i )=p 1>i X i +l3 2<i D i , (3{D i )=p Z!i X i +p i<i D i , I < i < n, (40) 

and (3(h) contains a monomial o,i 1 ,...,i n D % ^ , . . . ,DJ^ with a,i lt ...,i n ^ and i\ + 
. . . + i n = e, i.e., e = dcg Di D (3(h). After that one can find an automorphism 
7 such that j{X ) = X , 

j(Xi)= -Yi,i,3 X j, 7(A) = 7W D h l<i<n, (41) 

and (7 o (3)(h) contains a monomial aD e n with a coefficient / a £ F. Put 
a — 7 o p. We leave to prove the last assertion to the reader. The lemma is 
proved. 
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We apply the automorphism a. In what follows to simplify the notation we 
shall suppose without loss of generality that a = 1. So h contains a monomial 
aD e n with a coefficient O^oef, where e = deg ft.. It follows from here that 

deg Dji <*i,i = deg<5 M , deg Dn gjj = degg jd , h + 1 <j < k. (42) 

Let z = [z\, . . . , Zk) G A k be a solution of (|3D|) . Then (|4"2"f implies that one can 
uniquely represent 

z j= z 'j9j,j + X -.•- /) - ii + i<i<*. ( 43 ) 

0<s<deggj,j 

where ZpZj >a G h A, the degrees deg Dn Zj :S < for all Zi + 1 < j < k, < s < 
dc 

S-Di Again by (|42|) one can uniquely represent 

0<s<dcg 

where u",u' is G -A, the degrees deg^, tt^ s < for all 1 < i < I, < s < 
deg Di ^jj. Finally, by ([42l for all /i + 1 < j < fc, 1 < i < h, < r < deg Di gjj, 
one can uniquely represent 

0<r<dcg (5x,l 

where <5j, r ,i, $j,r,i,s G h A the degrees deg^ $j,r,i,s < for all considered j, r, z, s. 
Put 

I = { ti,r) : h + 1 < 3 < kkO < r < degg jd } , 
J = {(i,s) : 1 < i < h & 1 < s < deg <J M } . 

Therefore, 

Zi=~ Z 'jdj,i- X! Z 3,r8j,r,i + u", I < i < h, (44) 

ii + l<j"<fe U,r)ei 

z 3,rhr,i,» = u 'i,s, {t,s)&J. (45) 

0»ex 

Let us introduce new unknowns Zj <r , (j,r) G X. By (|43 j) (|43|) system ([37)) is 
reduced to the linear system 

Zj,T&i,T,i,8=u' ita , (i,s)ej. (46) 

(i,r)ez 

More precisely, any solution of system (|3"Tj) is given by (|4*3"|) . where G h A 
are arbitrary and Zj >r is a solution of system (|4"B1 over M. (we underline that 
here this solution Zj> may depend on D n although one can restrict oneself 
by solutions Zj >r which do not depend on D n ). Note that all 5j, r ,i,s and a 
are homogeneous elements of h A and there are integers d;, r , (j, r ) G 2, d^ s , 
(i, s) G J, p such that deg &j, r ,i,a = dj <r — d[ s and degw^ s = —d^ s + p for all 
(j, r) G I, {%, s) G J . This follows immediately from the described construction. 
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Now all the coefficients of system (|46l) do not depend on D n . As we have 
proved if the coefficients of (p?7| do not depend on X„ then the coefficients of 
(HU) also do not depend on X„, and hence in the last case they do not depend 
on X n , D n - 

If the coefficients of (|46j) depend on X n we perform an automorphism X n i— » 
D n D n i !■ — X n , X,- i— > Xj, D, I— ¥ Di, 1 < i < n — 1. Now the coefficients of 
system (|4"o]) do not depend on X n (but depend on D n ). After that we apply 
our construction recursively to system (|46p . 

The final step of the recursion is n = (although in the statement of theorem 
n > 1, see Section [T] we are interested only in Weyl algebras). In this case 
T = J = 0. Hence using (|44f for n = we get the required 2:* and J for n = 0. 

Thus, by the recursive assumption we get a particular solution Zj >r = z* r , 
(j, r) £ I, of system pp]l , an integer v\ (in place of v from assertion (A)) such 
that 

min {ordz* r } > min {ord u{ s } — v\ , (47) 

U,r)ex J ' (>,«)€J 

and a system of generators 

eI , 1 < a < /?, (48) 

of the module J' of solutions of the homogeneous system corresponding to (jUJ). 
Notice that if the coefficients of (|37|) do not depend on X n then J' is a module 
over the homogenization F[Xq,Xi, . . . , X n _i, £>i, . . . , D n _i] of the Weyl alge- 
bra of Xi, . . . , X n _i, Dt, . . . , D n _i. But obviously in the last case (pf8"|) gives 
also a system of generators of the M-module J" = h AJ' of solutions of the 
homogeneous system corresponding to (|4l)|) . Put 

Z* = - ^ z j,r$3,r,i + 11 i > 1 < * < il, 
(i,r)6I 

*i = E z i% £, n S «l + l<i<fc, 

0<s<deggj >: , 

* / * *\ 

Z — ^1 , • ■ ■ , Z fc J. 

Then z* is a particular solution of ([57|l . Put 

Za.,% = - y j z a ,j,r8j,r,i, l<i<h, 1 < ol < j3, 
(i,r)ez 

z aj = z a ^^ s D s n , h + l<j<k,l<a<{3, 

0<s<deggj t j 

z ^ h+jj i = 0, h + l<i,j <k, j % 

z f)-h+jj = 9j,j, h + l<j<k, 

Zf3-i 1+j A = -9j,%, l<i<h, h + 1 < j <k. 

Then J = Eko^+H! ^(*a,i, ■ • • , Hence ( k), 1 < a < 

[3 + k — li, is a system of generators of the module J. By |4T|) and the definitions 
of it', m" and s we have ordz* > ordu — vq — v\. Put v = vq + v\. 

LEMMA 10 All the degrees deg Sj^, deg gj^, deg8j^ t i, deg8j. r ^ tS and v, see 
above, are bounded from above by (nld)°^ , the degrees degtt^ are bounded from 
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above d' + (nld)°^ , the degrees degu", degu' is are bounded from above by 
d'(nld)°^. Further, all ord u'l , ordu^ s are bounded from below by ordu — v. 
Finally, in system |^6| ) the number of equations j^J is bounded from above by 
(nld)° {1) and the number of unknowns j^T is bounded from above by k(nld) 1 ■ ' . 

PROOF This follows immediately from the described construction. 

Let us return to the proof of Theorem [3] Applying Lemma [TU] and recursively 
assertions (A) and (B) for the formulas giving z* and J we get (A) and (B) 
from the theorem. The last assertion (related to the case when all bij and Uj 
do not depend on D n ) has been already proved. The theorem is proved. 

8 Proof of Theorem [J for Weyl algebra 

Let a be the matrix from Section [TJ We shall suppose without loss of generality 
that the vectors (a^i, . . . , fij,z), 1 < i < k, are linearly independent over the field 
F. We have degOij < d. This implies k < l{ d \l n )- 

Put the matrix b = h a. Let us define the graded submodulcs of h I 

Jo = h A{b 1A , b u ) + ... + h A(b kA , b k ,i), 
J v = J : (X») = {ze h A l : zX» e J }, v>\. 

We have the exact sequence of graded 'Vl-modules 

h A k -> J -> 0. 

Further, 3 V C J u+ i C h I for every v > and h I = \J U>0 J v - Since h A is 
Noetherian there is N > such that h I = Jm- So to construct a system of 
generators of I it is sufficient to compute the least N such that I = Jjv and 
to find a system of generators of Jn- 

LEMMA 11 h I = Jjy for some N bounded from above by {dl) 2 ° { ' . There is a 
system of generators b\, . . . , b s of the module J/v such that s and all the degrees 
deg&„, 1 < v < s, are bounded from above by (dl) 2 ( ' . 

PROOF Let us show that the module Jat+i C Jn for N >v. Let u G Jn+i- 
Consider system ([57]) . By assertion (A) of Theorem [2] there is a particular 
solution z* of (|37[) such that ordz* > 1. Hence u S XoJn C Jn- The required 
assertion is proved. Hence h I — J v . 

Let us replace in ((57]) (m, . . . , u{) by (UiXq, . . . , UiXq), where U\, . . . ,U\ are 
new unknowns. Then applying (B) from Theorem [2] to this new homogeneous 
linear system with respect to all unknowns U\, . . . ,Ui, Z\ , . . . , Z k we get the 
required estimations for the number of generators of J v and the degrees of these 
generators. The lemma is proved. 

COROLLARY 1 Let (of i 5 ■ ■ • ,a<i,i), 1 < i < I, be from the beginning of the 
section and the integer N be from Lemma Then for every integer m > the 
F-linear space 

A m+N (ai tl , . . . , aij) + . . . + A m+N (ak,i, a k,i) 3 I m - (49) 
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PROOF By Lemma E] we have (J ) m +N 3 X^(J N ) m = X^( h I) m . Taking 
the affine parts we get (pf§|). The corollary is proved. 

Now everything is ready for the proof of Theorem [T] By Lemma [11] and 
Lemma [T] there is a system of generators of the module gr(J) with degrees 
bounded from above by {dl) 2 ° { . By Lemma fT"2l from Appendix 1 the Hilbert 
function H(gt(I),m) is stable for m > (dl) 2 ° {n) . By flT0J) Section [2] the Hilbert 
function H(I,m) is stable for all m > (dl) 2 ° {n) . 

Consider the linear order < on the monomials from h A l which is induced by 
the linear order < on the monomials from A, see Section |H Then the monomial 
submodule C I C c A l is defined, see SectionEJ where C A = F[Xq, . . . , X n , D\, . . . , 
D n ] is the polynomial ring. By Section 2] the Hilbert function H( c I,m) is 

stable for all m > (dl) 2 ° l ' . Hence all the coefficients of the Hilbert polynomial 
of C I are bounded from above (dl) 2 ° ( ' . Therefore, according to (f5Tj) the module 
C I has a system of generators with degrees (d/) 2 " . This means, see Section HJ 
that the module Hdt( h I) has a system of generators with degrees (dl) 2 ° {n) . 
Therefore, the degrees of all the elements of the Janet basis of h I with respect 
to the induced linear order < are bounded from above by (dl) 2 ° {n) . Hence by 
Lemma [3] Section H] the same is true for the Janet basis of the module I with 
respect to the linear order < on the monomials from A 1 . Theorem \T\ is proved 
for Weyl algebra. 



9 The case of algebra of differential operators 

Denote by B = F(Xi, . . . , X n )[Di, . . . , D n ] the algebra of differential operators. 
Recall that A C B and hence relations (JT|) are satisfied. Further, each element 
/ £ B can be uniquely represented in the form 

/= E /, X f> l)J - 

jl,-,jn>0 3&1 

where all fj ly ...j n = fj G F(Xi, . . . , X n ) and F(Xi, . . . , X n ) is a field of rational 
functions over F. Let us replace everywhere in Section [1] and Section [2] A, 
XW 1 , degf = deg Xu Xn Du Dn f, dim F M, e Vtilj , f v>itj G F, (v,i,j), 

{i',f), by B, D\ degf = deg Di Dn /, <Km F (x 1 ,...,x n ) M > e v,j, fv,j 6 

F(Xi 7 . . . ,X n ), (v,j), j, j', j" respectively. Thus, we get the definition of the 
Janet basis and all other objects from Section Q] for the case of the algebra of 
differential operators. 

We define the homogenization h B of B similarly to h A, see Section[3l Namely, 
h B = F(X\, . . . , X n )[Xo, D\, . . . , D n ] given by the relations 

X..X, X..X,. D,l), I), I),, for all , 
DiXi - XiDi = X , 1 < i < n, A,/), /),.V, for all i^j. { } 

Further, the considerations are similar to the case of the Weyl algebra A with 
minor changes. We leave them to the reader. For example, Theorem [5] for 
the case of the algebra of differential operators is the same. One need only to 
replace everywhere in its statement A, h A and X n by B, h B and D n respectively. 
Thus, one can prove Theorem[T]for the case when A is an algebra of differential 
operators (but now it is B). Theorem [T] is proved completely. 
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One can consider more general algebra of differential operators. Let T be a 
field with n derivatives D\, . . . , D n . Then K n = !F[Dx, . . . , D n ) is the algebra of 
differential operators and similarly one can define its homogcnization h K n by 
means of adding the variable Xq satisfying the relations 

D%Dj = DjDi, X D i = D i X , A/ - /A = f Di X 

for all i,j and any element / G T where fu t G T denotes the result of the 
application of A to /. Following the proof of Theorem [1] one can deduce the 
following statement. 

REMARK 6 A similar bound to Theorem Q] holds for K n . 

Appendix 1: Degrees of generators of a graded 
module over a polynomial ring and its Hilbert 
function. 

We give a short proof of the following result, cf. pQ, [12], [H], 0]. 

LEMMA 12 Let I C A 1 be a graded submodule over the graded polynomial 
ring A — F[Xq, ■ ■ ■ , X n ], and I is given by a system of generators f±, . . . , f m of 
degrees less than d. Then the Hilbert function H(A /I,m) = dimpiA 1 / I) m is 
stable for m > {dl) 2 ° { + ' . Further, all the coefficients of the Hilbert polynomial 
of A 1 /I are bounded from above by (dl) 2 ° l + . 

PROOF Denote M = A 1 /I. Let L E F[X Q , . . . ,X n ] be a linear form in gen- 
eral position. Denote by K the kernel of the morphism M — > M of multiplication 
to L. We have K = {z G A 1 : Lz — J2 1<i<m fiZi, & z% G .4}. Hence solving 
a linear system over A, we get that K has a system of generators g\ , . . . , <? M 
with degrees bounded from above by (dl) 2 ' + . Let ?p be an arbitrary associ- 
ated prime ideal of the module M such that ^3 ^ (X , . . . , X n ). Since L is in 
general position we have L $ *p. Hence is not an associated prime ideal of 
K. Therefore, Km — for all sufficiently big N . So Xf gj G / for sufficiently 
big N and all i,j. Hence gj = Y,i<i< m Vj,ifi where y jti G F(Xi)[X , . . .,X n ]. 
Solving a linear system over the ring F(Xi)[Xo. . . . , X n ] we get an estimation 
for denominators from -F[^i] of all yj^. Since all gj and fi are homogeneous we 
can suppose without loss of generality that all the denominators are Xf . Thus, 
we get an upper bound for N. Namely, N is bounded from above by (d/) 20< ™ +1> . 
Therefore, the sequence 

-> M m -> M m+1 -> (M/LM) m+1 -> (51) 

is exact for m > (dl) 2 ° (n+1 ' . But M/LM = A 1 /{I + LA 1 ) is a module over 
a polynomial ring of F[X , ■ ■ ■ , X n ]/(L) ~ F[Xo, , . . , X n _i]. Hence by the 
inductive assumption the Hilbert function H(A l /(I + LA l ),m) is stable for 
to > (dl) 2 ° in) . Therefore, ([ST]) implies that the Hilbert function H(A l /I,m) is 
stable for m > {dl) 2 ° { + ' . 

Obviously for to < {dl) 2 °' + ' the values H(A l /I,m) are bounded from 
above by (gK) 2 ° ( ™ +1) . Hence by the Newton interpolation all the coefficients of 
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the Hilbert polynomial of A 1 /I are bounded from above by (dl) 2 ° ( + . The 
lemma is proved. 

We need also a conversion of Lemma [12] 

LEMMA 13 Let I C A 1 be a graded submodule over the graded polynomial 
ring A = F[Xq, . . . , X n ]. Assume that the Hilbert function H(A l /I,m) = 
dim p (A 1 /T) m is stable for m > D and all absolute values of the coefficients 
of the Hilbert polynomial of the module A 1 /I are bounded from above by D for 
some integer D > 1 . Then I has a system of generators f\ , . . . , f m with degrees 

PROOF Let us choose fx, . . . , f m to be the reduced Grobner basis of I with 
respect to an admissible linear order < on the monomials from A 1 , cf. the 
definitions from Section [1] and Section |U The degree of a monomial from A 1 is 
defined similarly to Section [1] and Section 2J We shall suppose additionally that 
the considered linear order is degree compatible, i.e., for any two monomials 
z\ , Z2 if deg z\ < deg zi then z\ < Z2 ■ For every z S A the greatest monomial 
Hdt(z) is defined. Further the monomial ideal Hdt(I) is generated by all Hdt(z), 
z G I. Now Hdt(/i), . . . ,Hdt(/ m ) is a minimal system of generators of Hdt(7) 
and deg fi — degHdt(/i) for every 1 < i < m. The values of Hilbert functions 
H(A l I 'Hdt(/),m) = H(A l /I,m) coincide for all m > 0. Thus, replacing / 
by Hdt(/) we shall assume in what follows in the proof that / is a monomial 
module. 

For every 1 < i < I denote by Ai C A 1 the z-th direct summand of A 1 . Put 
I L = I fl Ai, 1 < i < I. Then / ~ ®i<i<ili since / is a monomial module. 
Further, for every 1 < a < m there is 1 < i < I such that f a £ 2^. Let us 
identify Ai = A. Then Ii C A is a homogeneous monomial ideal. The case 
Ii = A is not excluded for some i. For the Hilbert functions we have 

H(A l /I,m) = H(A/Ii,m), m > 0. (52) 
i<i<i 

If (A/Ii)D = for some i then (A/Ii) m — for every m> D. In this case the 
ideal Ii is generated by ^Z Q<m<D {Ii)m- Hence in ([5*2")) for the values m > D one 
can omit this index i in the sum from the right part. Therefore, in this case the 
proof is reduced to a smaller I. So we shall assume without loss of generality 
that (A/h) D ¥=0,l<i<l 

Further, we use the exact description of the Hilbert function of a homoge- 
neous ideal, see [4] Section 7. Namely there are the unique integers fr^o > bi t \ > 
■ • • > bi t . n+ 2 = such that 

TT,,,-r \ fm + n+l\ yr-^ fm — bin + j — 1\ -_„. 

H(A/Ii,m) = ( - 1 - E ( 53 ) 

V 7 l<j<n+l 

for all sufficiently big m and 

6,.o = min{d : d > b hl &V m > d ^ holds}. (54) 

This description (without constants bifl) is originated from the classical paper 
|llj . The integers bi t o, . . . , bi :U +2 are called the Macaulay constants of the ideal 
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Jj. Besides that, 

fr(i,m) = H(A/Ii,m) - ( ^ J+1+ ]T ,lJ . J J>0 

V n ' l<j<n+l \ 3 ' 

(55) 

for every m > b^x, see [3] Section 7. By Lemma 7.2 [3] for all 1 < a < m if 
f a G Ii then deg/ a < &i.o- Hence it is sufficient to prove that all bi^, 1 < i < i, 
are bounded from above by Z) 20< + ' . 

By (l52l and (|55|) the coefficient at m™ -3 ', < j < n, of the Hilbert polyno- 
mial of A 1 1 1 is 



. .,. V 6i,„+l-j+ V V 7 — ry/i^ (56) 

v J ' l<i<l 0<v<j-l l<i<l y ' 

where ^ fij is an integer and fij V G Z[Z], < v < j — 1, is a polynomial 
with integer coefficients with deg fij_ v = j — v + 1. Moreover, and absolute 
values of all the coefficients of all the polynomials bounded from above 

by, say, 2°(" ). Denote bj = £ 

i<i<i — 3 ' — n + 2- By the condition of 
the lemma all the coefficients of the Hilbert polynomial of A /I are bounded 
from above by D. Hence from (|56p one can recursively estimate b n , . . . , b±. 



Namely, b n+1 -j = {2 n 'lD) 2 ° (j+1) , < j < n. Hence b x = {ID) 2 ^ L ' . Notice 
that bi t i < maxi<i<; 6^1 < b\ for every 1 < i < m. 

Now let m > maxi<i<; foj.i . By |53|) if /i(i, m) 7^ for some 1 < i < I then 
m < D, i.e., m is less than the bound D for the stabilization of the Hilbert 
function of A /I. Thus, fc^o < max{6j ) i, D} by (fS"4"|) . Hence 6^0 is bounded 
from above by (ID) 2 °' + ' . 

We have (>1/Ji)u + for every 1 < i < I. This implies H(A l /I,D) > I. 
Denote by Cj the j-th coefficient of the Hilbert polynomial of the module A / I. 
Now \cj\D j > l/(n + l) for at least one j. Hence D n+1 (n + 1) > I by the 

condition of the lemma. This implies that I 2 ( + ' is bounded from above by 
proved. 



D 2 . Therefore, bi t o is bounded from above by D 2 . The lemma is 



Appendix 2: Bound on the Grobner basis of a 
monomial module via the coefficients of its Hilbert 
polynomial 

Denote by C; = Z™ U • • ■ U Z™ the disjoint union of / copies of the semigrid 
Z™ = {(ii,...,i n ) : ij > 0, 1 < j < n}. A subset of C; which intersects 
each disjoint copy of I/L by a semigroup closed with respect to addition of 
elements from U] is called an ideal of C;. Any ideal I in Ci has a unique 
finite Grobner basis V = Vj, denote T = Ci \ I. Clearly, / corresponds to a 
monomial submodule in the free module (F[Xi, . . . ,X n ]) 1 . The degree of an 
element u = (fc;ii, ...,i n ) G C/, 1 < k < I is defined as \u\ — i± + ■ ■ ■ + in- 
The degree of a subset in C/ is defined as the maximum of the degrees of its 
elements. The Hilbert function Ht{z) equals to the number of vectors u G T 
such that |u| < z. Then Ht(z) = Eo<Km CsZS ' z — z o ^ or suitable zq, 
integers cq, . . . , c m where the degree m < n. Denote c = maxo< s < m |c s |s! + 1. 
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PROPOSITION 1 (cf. JBS, EM, W- The de 9 ree °f V does not exceed 
(cnY 

PROOF An s-cone we call a subset of a fc-th copy of Z™ in C; for a certain 
1 < k < I of the form 

P |.V„ /; Y, =l n - s } (57) 

for suitable 1 < j\, . . , ,j n - 8 < n. The degree of (1ST)) we define as \P\ = 
i± + • • • + i n -s (note that this definition is different from the one in [?]). By 
a predessesor of (|57p we mean each s-cone in the same fc-th copy of Z™ of the 
type 



— ii, . . . , Xj p l — ip—i, Xj p — i p 1, Xj p+1 — ip+i, . . . , Xj n _ B — *n-s} 

(58) 

for some 1 < p < n — s, provided that i p > 1. Fix an arbitrary linear order on 
s-cones compatible with the relation of predessesors. 

By inverse recursion on s we fill gradually T (as a union) by s-cones. For the 
base we start with s = m. Assume that a current union To C T of m-cones is 
already constructed (at the very beginning we put To = 0) and an m-cone of the 
form (fSTj) with s — m is the least one (with respect to the fixed linear order on 
m-cones) which is contained in T not being a subset of To. Observe that each 
predessesor of this m-cone was added to To at earlier steps of its construction. 
Since the total number of m-cones added to To does not exceed c m ml < c we 
deduce that the degree of every such m-cone is less than c rn ml (taking into 
account that the very first m-cone added to To has the degree 0). 

For the recursive step assume that the current To is a union of all possible 
m-cones, (m — l)-cones,...,(s + l)-cones and perhaps, some s-cones. This can be 
expressed as deg(Hr — Ht ) < s. Again as in the base take the least s-cone of 
the form ([5"7|) which is contained in T not being a subset of Tq. Observe that 
each predessesor of the type (|58p of this s-cone is contained in an appropriate 
r-cone Q, r > s, such that Q was added to To at earlier steps of its constructing 
and Q C {Xj p = i p — 1}. Hence 

\Q\>i p -l. (59) 

The described construction terminates when Tq = T. Denote by t s the 
number of s-cones added to Tq and by k s the maximum of their degrees. We 
have seen already that t m , k m < c. 

Now by inverse induction on s we prove that t Sl k s < {en) 2 ' ' . To this end 
we introduce a relevant semilattice on cones. Let C = {C aj ^} Qj /3, < (3 < 7 Q 
be a family of cones of the form (|5T[) where dimC Q ./5 = a. By an a-piece we 
call an a-cone being the intersection of a few cones from C. All the pieces 
constitute a semilattice C with respect to the intersection and with maximal 
elements from C. We treat C also as a partially ordered set with respect to the 
inclusion relation. Clearly, the depth of C is less than n. Our nearest purpose 
is to bound from above the size of C. For the sake of simplifying the bound we 
assume (and this will suffice for our goal in the sequel) that 7 Q < {en) 2 ( 
for s < a < m and 7 Q — when a < s, although one could write a bound in 
general in the same way. Besides that we assume that the constant in 0{. . .) is 
sufficiently big. In what follows all the constants in 0{. . .) coincide. 
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LEMMA 14 Under the assumption on the numbers j a < (cn) 2 , s < 

a < m of maximal elements of all dimensions from C, the number of a-pieces in 
L does not exceed (cn) 2 ( >+1 for s < a < m or (cn) 2 ° ( (s— w /j grl 
a < s. 

PROOF For each a-piece choose its arbitrary irredundant representation as 
the intersection of the cones from C. Let 6 be the minimal dimension among 
these cones. Then this intersection contains at most 8 — a + 1 cones. Therefore, 
the number of possible a-pieces does not exceed 

J2 M 2 ° (m " )( ^ a+1 \ 

max {a ,s} <~5 <rn 

that proves the lemma. 

Now we come back to estimating t Sl k s by inverse induction on s. Let in 
the described above construction the current To is the union of all added m- 
cones, (m — l)-cones,...,s-cones. Denote this family of cones by C and consider 
the corresponding scmilatticc C (see above). Our next purpose is to represent 
To as a Z-linear combination of the pieces from C by means of a kind of the 
inclusion-exclusion formula. We assign the coefficients of this combination by 
recursion in C. As a base we assign 1 to each maximal piece, so to the elements 
of C. As a recursive step, if for a certain piece P £ C the coefficients are already 
assigned to all the pieces greater than P, we assign to P the coefficient ep in 
such a way that the sum of the assigned coefficients to P and to all the greater 
pieces equals to 1. Therefore, we get 

T = ]T e P P 

Pec 

where the sum is understood in the sense of multisets. Hence 

H To (z)=^e P [ (60) 



Pec 



for large enough z. We recall that deg (Ht — Ht ) < s — 1. 

Now we majorate the coefficients \ep\ by induction in the semilattice C. The 
inductive hypothesis on t a < (cn) 2 ' 1 , s < a < m and Lemma fT4l imply that 

^2 l £p l - (cn) 2 ° <m " A) , s-l<X<m. 

dimP—X 

by inverse induction on A following the assigning ep. In fact, one could majorate 
in a similar way also X)dimP=A \ €p \ wnen A < s — 1, but we don't need it. The 
inductive hypothesis on k a < (cn) 2 ° l , s < a < m and (|60p entail that 
the coefficient of Ht (z) at the power z a does not exceed (cn) 2 ° ( \ s — 1 < 
a < rn (actually, due to the inequality deg(Hr — Ht ) < s — 1 the coefficients 
at the powers z a for s < a < m are less than c) . In particular, the coefficient at 
the power does not exceed (cn) 2 ° {m s+1> . Denote Ht — Ht = nz s ~ l + ■ ■ ■. 
By constructing To we add to it t s _i = n(s — 1)1 of (s — l)-cones, which justifies 
the inductive step for t s -i < (cn) 2 



2 0(m- S + l) 
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To conduct the inductive step for k s _i < (cn) 2 we observe that for 

each (s — l)-cone P added to T either every its predessesor is contained in a 
cone of dimension at least s, or some its predessesor is an (s — l)-cone as well. 
In the former case \P\ < (max s < Q < TO k a + l)(n — s + 1) (due to ([59]) ). while 
in the latter case |P| is greater by 1 than the degree of this predessesor, hence 
fc,s_i < (max s <„< m k a + l)(n — s + 1) + i s -i- Finally, exploit the inductive 
hypothesis for k m , . . . , k s , and the just obtained inequality on i s _i. 

To complete the proof of the proposition it suffices to notice that for any 
vector from the basis V treated as an 0-cone, each its predessesor of the type 
(155)) for s = is contained in an appropriate r-cone, whence the degree of V 
does not exceed (maxo< Q < TO k a + l)n again due to ([59")) (cf. above). 

Acknowledgement. The authors are grateful to the Max-Planck Institut 
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